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ABSTRACT 

We report on observations of the sixth accretion-powered milhsecond pulsar, ICR J00291-f 5934, 
with the Rossi X-Ray Timing Explorer. The source is a faint, recurrent X-ray transient initially 
identified by INTEGRAL. The 599 Hz (1.67 ms) pulsation had a fractional rms amplitude of 8% in 
the 2-20 keV range, and its shape was approximately sinusoidal. The pulses show an energy-dependent 
phase delay, with the 6~9 keV pulses arriving up to 85 fis earlier than those at lower energies. No 
X-ray bursts, dips, or eclipses were detected. The neutron star is in a circular 2.46 hr orbit with a 
very low-mass donor, most likely a brown dwarf. The binary parameters of the system are similar to 
those of the first known accreting millisecond pulsar, SAX J1808.4— 3658. Assuming that the mass 
transfer is driven by gravitational radiation and that the 2004 outburst fluence is typical, the 3-yr 
recurrence time implies a distance of at least 4 kpc. 

Subject headings: binaries: close — pulsars: individual (ICR J00291+5934) — stars: neutron — stars: 
low-mass, brown dwarfs — X-rays: binaries 
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1. INTRODUCTION 

The growing sample of accretion-powered millisecond 
X-ray pulsars divides naturally into two groups. Three 
of the five sources known to date (XTE J175 1— 305. 
[Mark wardt et al. 2002; XTE J0929-314, Galloway et all 
I200I and XTE J1807-294, Markwar dt et all l2003() 
are in ultracompact binaries with orbital periods of 
« 40 min. The Roche lobes in such tiny bina- 
ries cannot contain a main-sequence companion, in- 
dicating that the mass donors are highly evolved 
and H -poor l^ clson et al. 198^ iNelson fc RappaportI 
l2003t iDelove fc Bildsten .20031). The othe r two 
(SAX J1808.4-3658, Wiinands fc van der KM 11993 



[Chakrabartv fc Morgan' '199^, and XTE J1814-338, 
iMa rkwardt & Swank 2003) have orbital periods of 2.01 
and 4.28 hr, respectively, and have H-rich donors (likely 
a brown dwarf in the case of SAX J1808.4— 3658; 
iBildsten fc Chakrabartvll2001[l . The latter two sources 
have also both exhibited thermonuclear (type-I) X-ray 
bursts, whereas the three ultracompact binaries have not. 
All five pulsars are soft X-ray transients, with outburst 
durations of order weeks and recurrence times of order 
years. It remains unclear why millisecond pulsations are 
easily detectable in these five sources, but not in the over 
50 other neutron star s in low- mass X -ray binaries (e.g. 
iCumming et al.ll2001l: FRtarchuk et a l. "2002) . 

The 2002 October launch of th e INTEGRAL gamma- 
ray mission jWinkler et al.ll2003|) . with its wide field of 
view and good sensitivity to hard X-ray sources, has 
provided a new avenue for detections of X-ray tran- 
sients. INTEGRAL discovered the new X-ray transient 
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ICR J00291+5934 (I = 120?1, b = -3?2) on 2004 Dec. 
2 (MJD 53341; Eckert et al. 2004). FoUowup Rossi X- 
ray Timing Explorer (RXTE) observat ions revealed pul- 
sat ions at a frequency of 598.88 Hz ijMarkwardt et af] 
I200 4b1. and exhibiting a sinusoidal f requency modula- 
tion indicative of a 147.4 min orbit fMark wardt et af] 
[^04a). The detection of pulsations motivated extensive 
multiwavelength foUowup. An R w 17.4 candidate op- 
tical counterpart was identified within the INTEGRAL 
error circle, at R.A. = 00*'29™03';06, deck = -H59°34'19"0 
(02000.0, la undertainty radius 0'.'25; iFox fc Kulkml 
|2004|) . Optical spectroscopy of the candidate revealed 
weak Hell and Ha emis sion, supporting its association 
with ICR J0029H-5934 l|R,oelofs et al.ll200l . Radio ob- 
servations also revealed evidence for vari able ernission 
consistent with the counterpart position ijPoolevI l2004t 
Fender et al. 2004). 

In this letter we describe the detailed analysis of the 
RXTE observations following the discovery. 

2. OBSERVATIONS 

We obtained a series of pointed RXTE observations of 
ICR J0029H-5934 between 2004 December 3-16 (MJD 
53342-53355), with a total exposure of approximately 
295 ks. Our analysis is primarily based on data from the 
RXTE Proportional Counter Array (PCA lJahoda et alJ 
Il99fi), which consists of five identical gas-filled propor- 
tional counter units (PCUs) sensitive to X-ray pho- 
tons in the 2.5-60 keV range, with a total effective 
area of « 6000 cm^. The data were collected either 
in GoodXenon or generic Event modes (in addition to 
the standard data modes). GoodXenon records the ar- 
rival time (1 IIS resolution) and energy (256 channel res- 
olution) of every unrejected photon, while the chosen 
Event configuration has 125 /is time resolution and 64 
energy channels. We analysed FITS production data us- 
ing LHEASOFT version 5.3.1 (2004 May 20) and esti- 
mated the background flux using the "CM" bright-source 
model for PGA gain epoch 5 (from 2000 May 13). The 
photon arrival times at the spacecraft were converted 
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Fig. 1.— Pulse arrival time delays in IGR J00291+5934 due to 
the 2.46 hr binary orbit, with respect to a constant pulse frequency 
model. The measured time delays are plotted as filled circles, and 
the symbols with error bars (representing the la uncertainties) 
show the model residuals (on a 100 X expanded scale, right-hand 
j/-axis) . 

to barycentric dynamical times (TDB) at the solar sys- 
tem barycenter using the Jet Propulsi on Laborato ry DE- 
200 solar system ephemeris (Standis h et al.lll992i) along 
with a spacecraft ephemeris and RXTE fine clock cor- 
rections. The positio n adopted was that of the proposed 
optical counterpart fF ox fc KulkarnI '2004*1. Data from 
the 20-200 keV High-E nergy X-ray Timing Experiment 
(HEXTE; Gruber et al HTooq) were also used to charac- 
terize the hard X-ray spectrum. 

3. ANALYSIS AND RESULTS 

For our timing analysis, we selected photons from the 
top layer of each PCU in the energy range 3-13 keV (ab- 
solute channels 7-30) in order to maximize the pulsed 
signal to noise ratio. We binned the arrival times for 
these photons into 2~^^ s («0.12 ms) samples. A 599 Hz 
pulsed signal was easily detectable in these data, with a 
fractional rms amplitude of between 7-9%. The rms am- 
plitude decreased throughout the outburst, gradually to 
6% on December 11 (MJD 53350) and then more rapidly 
down to a minimum of « 2% by December 14. The pul- 
sations became undetectable (< 1% rms) from December 
15 (MJD 53354) onwards, and the residual flux which was 
detected after this date likely arose from other sources in 
the RXTE field of view. The pulse frequency, measured 
every 64 s, exhibited a sinusoidal modulation attributed 
to Doppler shifts from the binary orbit. Wc fitted a pre- 
liminary orbital model to the frequency measurements 
between December 3-6 (MJD 53342-53345), and then 
used this prel iminary model in a more pre cise pulse phase 
analysis (e.g. iManchester & Tavlojll977!) . We first cor- 
rected the time series using the preliminary orbit. We 
then folded 64-s intervals of data on the pulsar period 
estimated from the frequency fit, and determined the 
epoch of maximum flux within each interval. We fitted 
the residuals (measured pulse arrival time minus the pre- 
dicted arrival time) to compute diffe rential corrections t o 
our Keplerian orbit model (see, e.g. lDeeter et alJll98lD . 
and repeated the process with the corrected model pa- 
rameters until no further reduction in the rms residual 



error was achieved. No significant eccentricity was de- 
tectable. Our best-fit orbit and spin parameters for the 
pulsar are given in Table 1. The pulse time delays due 
to the orbit are shown in Fig. ^ 

We extracted light curves in 12 energy bands between 
2-20 keV with 0.12 ms time bins from the observations 
between December 3-6 (MJD 53342-53345). We then 
corrected these light curves for orbital motion based on 
the best-fit orbital solution, and folded on the pulsar pe- 
riod to measure the energy dependence of the pulsations. 
The fractional rms pulsed amplitude decreased with en- 
ergy, from 10% in the 2-3 keV range to 6.5% between 
7-16 keV, dropping further to 5% above 16 keV (Fig. 
121 upper panel). The profile was approximately sinu- 
soidal at all energies. Of the higher harmonics, only the 
second harmonic (1198 Hz) was detected at > 3cr con- 
fidence in more than one energy band (although not in 
all bands), at fractional rms amplitudes of a few tenths 
of a percent. Upper limits on the fractional amplitudes 
for the third and higher harmonics were of similar mag- 
nitude. The pulse phase had an energy dependence. The 
hard photons arrived earlier than the soft, with a maxi- 
mum lead of around 90 /zs (equivalent to 5% of the pulse 
period) at 6-9 keV (Fig. 13 lower panel). In contrast 
to other accretion-powered millisecond pulsa rs, where 
the hard lead satura tes above 10 keV ijCui et al. 199^ 
iGallowav etUI I2002D . in IGR J0029H-5934 we found 
evidence for a slight decrease in the lead above 9 kcV. 
We found no evidence for modulation of the X-ray in- 
tensity at the orbital period, ruling out eclipses or dip- 
ping behavior. We searched for kilohertz quasi-periodic 
oscillations (kHz QPOs) in the range 50-4000 Hz, aver- 
aging FFTs from 128 and 512-s segments of data over 
entire observations and entire days, but found no sig- 
nificant peaks. We estimate a typical 3(t upper limit 
for features with FWHM of up to 100 Hz of < 1% 
rms (compare to the 5 -10% rms kHz QPOs measured 
in SAX J1808.4-3658: IWiinands et af]l200l . We also 
searched each observation of IGR J00291-I-5934 for ther- 
monuclear X-ray bursts, but found none^. 

We fitted HEXTE and PGA spectra, the latter ex- 
tracted separately for each PGU (excluding PGU #0, 
which has lost its propane veto layer) between 2.5- 
200 keV. We obtained a good fit with a model consisting 
of a power law with spectral index 1.7-1.8, attenuated by 
neutral absorption with a column depth of ~ 10^^ cm~^. 
This is only sligh tly higher than the surv ey value of 
4.7 X 10^^ cm-2 (iDickev fc LockmanllT99(H) . A Gomp- 
tonization model ljTitarchumT994() gave a marginally 
better fit to the highest signal-to-noise spectra but the 
scattering electron temperature (which gives rise to a 
high-energy roll off in the spectra) could not be mea- 
sured, save that it was > 60 keV. This implies an opti- 
cal depth for scattering (assuming spherical geometry) of 
T — 0.8. We also found evidence for neutral Fe Ka line 

^ A burst was detected by both PGA and HEXTE at 18:34:21 UT 
during an observation of IGR J00291-I-5934 on 2004 December 
12. However, this event had a very unusual energy spectrum and 
was probably a bright gamma-ray burst observed through the side 
of the RXTE detectors. The burst was also detected by INTE- 
GRAL ( SPI/AGS monitor trigger 4779) and the Mars Odyssey 
mission |HurlfiJ!_et al."2002'), and the derived interplanetary net- 
work (IPN3; see \protect ittp://www.ssl. berkeley.edu/ipn3/ I po- 
sition annulus excludes IGR J00291+5934 as the origin (K. Hurley, 
pers. comm.). We therefore designate this event as GRB041212. 
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Fig. 2. — Top panel Fractional rms amplitude of pulsations in 
IGR J00291+5934 as a function of photon energy. Bottom panel 
Pulse arrival time as a function of energy, measured relative to the 
pulsation in the 2-3 keV energy band (negative values correspond 
to earlier arrival times). In both panels, vertical error bars indicate 
the la uncertainties 

eraission at 6.4 keV, with equivalent width « 50 eV. 
While the line emission could potentially originate from 
other sources within the 1° RXTE field of view, or from 
diffuse Galactic emission, a Chandra observation on De- 
cember 14 also found evidence for Fe line emission in the 
zeroth order spectrum (M. Nowak, pers. comm.). In 
order to minimise the fit residuals for the PCA data we 
also added two edge features, one at 7.12 keV (Fe neutral 
K-edge) and one at 4.7 keV (instrumental xenon edge). 

An X-ray intensity history of the 2004 outburst is 
available both from the PCA observations and from 
the RXTE All Sky Monitor (ASM; Lcvinc ct al. 1996); 
these are shown in Fig. O The flux measured by 
the PCA decreased from 1.14 x 10~^ ergscm~^s~^ to 
< 10"^° ergscm-^s^i (2.5-25 keV) between 2004 De- 
cember 3-15 (MJD 53342-53354). During the dis- 
covery observations by INTEGRAL on December 2 
(MJD 53341) the JEM-X instrument measured a 3- 
10 keV flux of 23 ± 5 mCrab, while the ISGRI instru- 
ment measured 55 ± 5 mCrab in the 20-60 keV band 
(jEckert et alJl200^ . For an absorbed power law with 
riH = 10^^ cm~^ and photon index F = 1.7 (where 
dN/dE oc E^^), this translates to a 2.5-25 keV flux 
of 1.0-1.1 X 10"^ ergscm-^s"^. Thus, the peak PCA 
flux measurement of (1.14 ± 0.06) x 10~^ ergscm^^s^^ 
(2.5-25 keV) on December 3 (MJD 53342) is hkely a 
good measure of the peak outburst flux. We calculated 
a correction factor of 2.54 from the broadband spectral 
fits to apply to the 2.5-25 keV fluxes, in order to es- 
timate the bolometric flux (here taken as the flux from 
0.1-200 keV). Assuming the outburst commenced on De- 



cember 2 (MJD 53341), we estimated a total fluence for 
the outburst of 1.8 x 10~^ ergscm"^. 

The RXTE mission-long (1996-2004) 2-10 keV flux 
history of IGR J00291-I-5934 produced from ASM ob- 
servations confirmed the current outburst, as well as 
revealing evidence for two previous instances of activ- 
ity, on 1998 Novemb er 26-28 and 2001 September 11- 
21 ijRemillardI l2004|) . These results indicate that the 
transient appears regularly, roughly every three years, 
similar to the 401 Hz pulsar SAX J1808.4-3658, which 
exhibits outbursts every w 2 yr (e.g. Wlina nds 2003). 
Combined with the estimated 3-yr recurrence time of 
the outbursts in IGR J00291-I-5934, and assuming simi- 
lar fluences, we estimate a long-term time-averaged flux 
for IGR J00291-t-5934 of 1.8 x 10^" ergscm^^s-i. 

4. DISCUSSION 

With a spin frequency of 599 Hz, IGR J00291-f5934 
is the fastest known accretion-powered millisecond pul- 
sar. Although it is neither the fastest known pul- 
sar (cf. radio pu lsar PSR B1937+21, v = 641 Hz; 
IBacker et alllTgp) nor the fastest spinning neutron star 
in an X-ray binary (cf. nuclear-powered millisecond pul- 
sar 4U 1608-52, V = 619 Hz; Hartman et al. 2005, 
in prep.), it is part of a growing population of pulsars 
found in the 600-650 Hz range. These objects under- 
score the statistically signiflcant absence of pulsars with 
spins faster than 6 50 Hz which has been noted p revi- 
ouslv l|Chak rabart v"eraIll2003t IChakrabartvll2005j) . It 
remains unclear whether this cutoff in the spin distri- 
bution reflects the dist ribution of neutron s tar mag- 
netic fleld strengths fe.g. lWhite & Zh ang 1997) or angu- 
lar momentum losses fr om gr avitational radiation (e.g. 
IWagonei<IT98l IBildstenl(T99? 
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Fig. 3.— X-ray intensity of IGR J00291-I-5934 throughout the 
2004 outburst. The 1-day averaged ASM intensity values and their 
Icr confidence intervals are shown as shaded regions (left-hand y- 
axis) while the integrated 2.5-25 keV PCA flux measurements are 
overplotted as open squares (right-hand j/-axis). The Icr uncertain- 
ties on the PCA flux measurements are shown, but are generally 
smaller than the symbols. The estimated 2.5-25 keV flux during 
the intitial INTEGRAL detection is indicated by the open dia- 
mond. 
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In many respects, IGR J00291+5934 closely resembles 
the best-known millisecond pulsar, SAX J1808.4— 3658. 
The mass function in IGR J00291-I-5934 implies a min- 
imum companion mass of Mc = O.O39M0 (for a 1.4Af0 
neutron star and binary inclination i = 90°, although 
the lack of X-ray eclipses or dips implies i < 85°). For 
an isotropic a priori distribution of inclinations, we find 
an upper limit on the companion mass of O.IGMq (95% 
confidence, assuming = 2Mq). The inferred mass- 
radius relation intersects that of a hydrogen main se- 
quence star at a companion mass of O.24M0 (Fig. 21, 
which requires an improbably small inclination of 10°. 
This indicates that, hke SAX J1808.4-3658, the mass 
donor in IGR J00291-I-5934 is a hot brown dwarf, likely 
heated by low-level X-ray em ission during quiescence 
l|Bildsten fc Chakrabartvl200H) . For very low-mass com- 
panions, the mass transfer rate driven by gravitational 
radiation is 



AfcR ^ 3.3x10" 



0.039M( 



1.4M, 



2/3 



Pn 



2.46 hr 

(1) 

where Mc is the mass of the donor and the mass of 
the neutron star. Using the long-term averaged fiux of 
1.8 X 10""'^^ ergscm~^s~^, we derive a minimum distance 
of roughly 4 kpc (for a companion mass of 0.039Mq). 
This would place the system at least 200 pc away from 
the Galactic plane, likely within the disk. We note that 
the Mgr above refers to the mass transfer rate averaged 
over a much longer time scale than the RXTE mission 
lifetime, so that this distance estimate must be viewed 
with some caution. 

The peak inferred accretion rate for a distance of 4 kpc 
is > 0.04d|j,p^ MEdd (where MEdd - 2 x IQ-^ Mq yr"! 
is the hydrogen Eddington limit for spherical accre- 
tion onto a 1.4M0 neutron star). This is comparable 
to the maximum reached by SAX J1808.4— 3658 dur- 
ing its 2002 October outburst, when four thermonu- 
clear (type-I) X-ra y bursts were detected by RXTE 
l|Ghakrabartv et all l2003l) . The minimum burst in- 
terval in SAX J1808.4-3658 was 21 hr. If the dis- 
tance to IGR J0029H-5934 is greater than 4 kpc, the 
accretion rate would also be larger by a factor of 
((j/4 kpc)^, and we would e xpect more frequent bursts 
(e.g. iNaravan fc Hevll 12003(1 . In this case, however, we 
would expect to have detected at least one burst in our 
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Fig. 4.— Possible companion types for IGR J00291-I-5934. The 
dark solid curve shows the allowed mass-radius relation based on 
the RXTE observtions. The other curves show theoretical mass- 
radius relations for cold He white dwarfs {dotted curve), low- mass 
hydrogen main sequence stars {dashed curve), and brown dwarfs 
{dot-dashed curves) of ages 0.1, 0.5, 1.0 and 5.0 Gyr from top 
to bottom. The young (0.1 Gyr) brown dwarf model should be 
roughly representative of a bloated, X-ray-heated brown dwarf 
donor (see Bildstcn & Chakrabartv 2001) 

RXTE observations, since the observational duty cycles 
during the 2002 October outburst of SAX J1808.4-3658 
and the 2004 December outburst of IGR J00291-h5934 
were comparable. The non-detection of thermonuclear 
bursts from IGR J00291-h5934 thus suggests that the 
source distance is not significantly greater than 4 kpc. 
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TABLE 1 

Orbit and Spin Parameters for IGR J00291+5934 



Parameter 



Value 



Orbital period, Po^b (s) 

Projected semimajor axis, Oxsini (light-ms) 
Epoch of 90° mean longitude, T^/a (MJD/TDB) 
Spin frequency, vq (Hz) 
Mass function fx (^©) 
Eccentricity, e 

Spin frequency derivative, \u\ (Hz s~^) 



8844.092(6) 
64.993(2) 
53345.1875164(4) 
598.89213064(1) 
2.81311(7) X 10~5 
< 2 X 10-* (3(t) 
< 8 X 10-" (3(j) 



xVdof 



980/748 



Note. — These paramet ers assume the position of the opti- 
cal counterpart identified bv lFox fc K ulkarni (2004). Our timing 
analysis is based on the first 50 ks of data, between 2004 Decem- 
ber 3-6 (MJD 53342-53345). 
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